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Abstract 

An approximate procedure for performing nonperturbative calculations in quantum field theories is pre- 
sented. The focus will be quantum non-Abelian gauge theories with the goal of understanding some of 
the open questions of these theories such as the confinement phenomenon and glueballs. One aspect of 
this nonperturbative method is the breaking down of the non-Abelian gauge group into smaller pieces. For 
example 5(7(2) —*U(1) + coset or SU(3) — > SU(2) + coset. The procedure also uses some aspects of an old 
method by Heisenberg to calculate the n-point Green's function of a strongly interacting, non-linear theory. 
Using these ideas we will give approximate calculations of the 2 and 4-points Green's function of the theories 
considered. 



1 Introduction 

Modern perturbative quantum field theories (QFT) are based on the concept of a particle or field quanta. This 
perturbative paradigm has given rise to the most accurate description of some aspects of the physical world. 
The standard example being the theoretical calculation of the anomalous magnetic moment of the electron 
where theory and experiment are in agreement to the order of 10 -9 . Mathematically the elementary particles 
that arise in perturbative quantum field theory are harmonic excitations of fundamental fields. The quanta 
are defined through the creation and annihilation operators, a* and a, of the "second-quantized" theory. For 
instance, in quantum electrodynamics (QED), the entire electromagnetic field can be seen as a collection of 
superposed quanta, each with an energy Lu k ■ The Hamiltonian operator of the electromagnetic field (omitting 
the zero-point energy) can be written 

H = J2^ k , (1.1) 

k 

where 

N k = a\a k , (1.2) 

is the "number operator", i.e. giving the number of quanta with a specific four-momentum k when operating 
on a free state 

N k \—n k ...) = n k \...n k ...). (1.3) 
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where nk is a positive integer, the number of quanta with that particular momentum. The energy in the 
electromagnetic field is thus the eigenvalue of the Hamiltonian (1). The reasoning for fermion fields is the same, 
but then the number of quanta in any given state can be only or 1 due to Fermi statistics. 

Usually, in quantum field theories (for example, in QED) "fields" and "particles" are used interchangeably. 
This is because the almost universal usage of Feynman diagrams gives the false impression that a field is always 
equivalent to a "cloud" of quanta. The Feynman diagram technique (based on propagators and vertices) can 
be justified as an approximation for mildly nonlinear theories (i.e. in the weak coupling limit) but breaks down 
for strongly coupled Abelian and non-Abelian theories. 

The physical reason of the appearance of quanta is that in linear theories (for example, in electrodynamics) 
a general solution can be obtained by making the Fourier expansion: 

r 3 

A unear = / ^ k £ [ 0fc (A) E/J (*, A) e -*« + a\ (A)e* (*, A)e**] , (1.4) 

J A=0 

where e M is the polarization vector. 

However, for a theory based on a non-Abelian group, like quantum chromodynamics (QCD), this can no 
longer be done in the strong coupling limit, due to the nonlinear nature of Yang-Mills equations 

A bQCD^ I d ^ k ^[al(X)e^kA)e- ik -- + al\X)e;(k,X)e ik -% (1.5) 

Formally one can change this inequality to an equality but this is only a good approximation in the limit when 
the coupling is weak - g s <C 1. Due to the nonlinear character of non-Abelian gauge theories the superposition 
principle does not hold in general. There are plane wave solutions in non-Abelian gauge theories [1] however 
one can not superpose different plane waves as is done in (|1.4fl . Thus, color fields can be represented by 
harmonic oscillators (gluons) only in the limit when the strong interaction coupling constant tends to zero, 
g s — ► (or equivalently when Q 2 — ► oo because of asymptotic freedom). There are certain nonperturbative 
field configurations (for example, the hypothesized color electric flux tube that is thought to stretch between 
quarks and lead to the dual superconducting picture of confinement) in which the field distribution can not be 
explained as a cloud of quanta. The fields are there, but their quanta - gluons and quarks - are relevant only 
when probed at sufficiently short distances (or large momentum) . One way of looking at this situation is that 
the fields of these nonperturbative configurations are split into ordered fields (for example, the fields inside flux 
tube stretched between quark and antiquark) and disordered fields (the fields outside the flux tube). These 
field can not be interpreted as a (perturbative) cloud of quanta. In such situations the fields play a primary 
role over the particles 

The need for nonperturbative techniques in strongly interacting, nonlinear quantum field theories is an old 
problem that has been around since the beginning of the study of quantum fields. Much effort has gone into 
trying to resolve this puzzle. The different approaches that have been tried include: (i) lattice QCD [3]-JE|i (ii) 
the dual Meissner effect in the QCD- vacuum (hi) instantons [221 |23j . (iv) path integration |24 | -|29 j . 

(v) analytic calculations [3DI ? (vi) Dyson-Schwinger equations [31] (vii) and various other methods P2" ] -[4"? | . 
Despite this the problem is not yet fully resolved. In the 1950's Heisenberg |33| studied a nonlinear spinor 
field and worked out nonperturbative techniques for quantizing the nonlinear spinor field. In this method 
one writes down an infinite system of equations which connects together all the n-point Green's functions of 
the theory (this can be compared to the infinite number of Feynman diagrams which must, in principle, be 
calculated for a given process in pertubative quantum field theory) . In order to solve this system of equations 
one must find some physically reasonable approximations for cutting off this system of equations at some point 
so that one reduces the infinite system of equations into a finite system. Some preliminary work |5(J)-|53| has 
been done in applying these nonperturbative quantization techniques, in conjunction with ideas from Abelian 
projection, to study quantum SU(2) and SU(3) gauge theories, and also to examine the possible application of 
these techniques to the high-T c superconductors of condensed matter physics. The advantage of this approach 
is that it is nonperturbative. The problem is to find physically reasonable motivations for cutting off the infinite 
equations set for Green's functions. 

The preliminary work in this nonperturbative approach was devoted to a group of classical singular solutions 
of the Yang-Mills field equations, and their quantization [SD]-[S3j. In these papers some singular classical 
solutions were found, and they were then quantized using a nonperturbative method with the approximation 
that some of the degrees of freedom of the non-Abelian gauge theory could be treated as almost classical 
degrees of freedom (this was the approximation which cut off the infinite system of equations). After the 
quantization procedure was applied the bad asymptotic behavior (the divergence of some components of the 
Yang-Mills fields) of the original classical solution was modified into acceptable asymptotic behavior. Physically 
the nonperturbative quantization smoothed out the bad classical behavior of the solutions. 
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In later work [S^ the role of the other, non-classical degrees of freedom was investigated. In Ref. [SI] the 
above approach was applied to the SU(2) non-Abelian gauge theory. The assumption was made that the SU(2) 
gauge field could be split into components which had different qualitative behavior: one component, A^ , 
was taken to be in an ordered phase, while the remaining two components, A 1 ^ 2 , were in a disordered phase. 
Physically the ordered phase described the U(l) electric field in a hypothesized flux tube stretched between 
sources (i.e quarks) of this field, and the disordered phase described a condensate which "pushes out" (i.e. 
exhibits the Meissner effect) the ordered phase. The calculations of this [HI] indicated that such a scenario did 
occur. In this paper the dynamics of the disordered phase were "frozen" or taken to be nondynamical. 

In the present paper we would like to bring together and extend several results of these analytical calculations 
in nonperturbative QFT. Two general types of calculations will be presented: (i) the decomposition of the SU(3) 
gauge field components into SU(2) + coset components. This will lead to flux tube solutions which are important 
in the dual superconducting model of confinement, (ii) The simplification of the SU(3) quantized fields using a 
scalar field. This will lead to a simple model for a glueball. 

2 Heisenberg quantization method 

2.1 Initial Heisenberg idea for the quantization 
of a nonlinear spinor field 

Heisenberg's basic idea is that the n-point Green's functions can be found from some infinite set of coupled 
differential equations, which are derived from the field equations for the field operators. As an example we show 
how this method of quantization works for a spinor field with a nonlinear self interaction 0S| 0§| ■ The basic 
equation for the spinor field is : 

-fd^x) - l 2 Q[i>(x)(^(x)tP(x))} = (2.1) 

where 7 M are Dirac matrices; ip(x),ip(x) are the operators of the spinor field and its adjoint respectively; 
Sstyfipip)] = ip(ipip) or ■07 5 ('07 5 ^) or ipj^^^/^ip) or 

The constant I has units of length, and sets the scale for the strength of the interaction. 
Heisenberg emphasized that the 2-point Green's function, G2(x2,xi), in this theory differs strongly from the 
propagator in a linear theory in its behavior on the light cone : in the nonlinear theory G 2 (x 2 ,x\) oscillates 
strongly on the light cone in contrast to the propagator of the linear theory which has a 5-like singularity. 
Heisenberg defined r functions 

T(x lX2 ...\ yi y 2 ...) = (0\T$( Xl )ij>(x 2 ).4(yS(y2)---]\®) (2-2) 

where T is the time ordering operator; |$) is a state for the system described by Eq. (|2.1() . Eq. (|2.2() establishes 
a one-to-one correspondence between the system state, |$), and the set of functions r. This state can be defined 
using the infinite function set of Eq. (|2.2(l . Using equation 12.1J1 and 12. 2J) we obtain the following infinite system 
of equations for the r's 

d 

l~ 2 l( r \-^-T(xi...x n \y 1 ...y n ) = ^s[T(x 1 ...x n x r \y 1 ...y n y r )} + 

Ox (r) 

5(x r - y 1 )T(xi...x r -xx r+1 ...x n \y 2 ...yr-iyr+i—yn) + 

5(x r - y 2 )T(xi...x r -ix r+ i...x n \yiy2...yr-iyr+i--yn) + ■■• (2.3) 

Eq. 1)2.311 represents one of an infinite set of coupled equations which relate various orders (given by the index 
n) of the r functions to one another. To make some head way toward solving this infinite set of equations 
Heisenberg employed the Tamm-Dankoff method whereby he only considered r functions up to a certain order. 
This effectively turned the infinite set of coupled equations into a finite set of coupled equations. 

Heisenberg used the procedure sketched above to study the Dirac equation with a nonlinear coupling. Here 
we apply this procedure to nonlinear, bosonic field theories such as QCD in the low energy limit. In particular 
we will apply this method to several solutions of the Yang-Mills field equations which have divergent fields 
and infinite energy. By making certain assumptions analogous to the Tamm-Dankoff cut-off, the unphysical 
asymptotic behavior of these classical Yang-Mills solutions will be "smoothed" out. 

2.2 Heisenberg quantization for QCD 

In this subsection we will apply a version of Heisenberg's quantization method sketched above to QCD. The 
classical SU(3) Yang-Mills equations are 

d v T BiiV = (2.4) 
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where T B V = d^A B - d v A B + gf BCD A°A° is the field strength; B, C, D = 1, . . . , 8 are the SU(3) color indices; 
g is the coupling constant; f BCD are the structure constants for the SU(3) gauge group. In quantizing the 
system given in Eqs. Ij2.4[l - via Heisenberg's method one first replaces the classical fields by field operators 
A B — * A B . This yields the following differential equations for the operators 

d v T BliV = 0. (2.5) 

These nonlinear equations for the field operators of the nonlinear quantum fields can be used to determine 
expectation values for the field operators A B , where (■ • • ) = (Q\---\Q) and \Q) is some quantum state). 
One can also use these equations to determine the expectation values of operators that are built up from the 
fundamental operators A B . For example, the "electric" field operator, £f = d Af — d z A B + gf BCD A$ A® 

giving the expectation (Eg). The simple gauge field expectation values, (A^x)), are obtained by average Eq. 
(|2.5|l over some quantum state \Q) 



Q 



d v T B ^ 



Q)=0. (2.6) 



One problem in using these equations to obtain expectation values like (A B ), is that these equations involve 
not only powers or derivatives of (A B ) (i.e. terms like d a (A B ) or d a dp(A B )) but also contain terms like 
Gffi = (A B A°). Starting with Eq. one can generate an operator differential equation for the product 

A B A„ thus allowing the determination of the Green's function Q B „ 

A B (x)d yv T B ^(x) Q) = 0. (2.7) 



Q 

However this equation will in turn contain other, higher order Green's functions |5l)j -|52). Repeating these steps 
leads to an infinite set of equations connecting Green's functions of ever increasing order. This construction, 
leading to an infinite set of coupled, differential equations, does not have an exact, analytical solution and so 
must be handled using some approximation. 

Operators which are involved in equation (|2.5|l are only well determined if there is a Hilbert space of 
quantum states. Thus we need to look into the question of the definition of the quantum states \Q) in the above 
construction. The resolution to this problem is as follows: There is an one-to-one correspondence between a 
given quantum state \Q) and the infinite set of quantum expectation values over any product of field operators, 
G™::(xx,x 2 ...) = (Q\A™(xi)A™(x 2 ) ...\Q). So if all the Green's functions - £r™".7(xi, z 2 ■ ■ ■) - are known then 
the quantum states, \Q) are known, i.e. the action of \Q) on any product of field operators AV? (xi)A™(x2) ... is 
known. The Green's functions are determined from the above, infinite set of equations (following Heisenberg's 
idea). 

Another problem associated with products of field operators like A™(x)A"(x) which occur in Eq. H2.5JI is 
that the two operators occur at the same point. For non-interacting field it is well known that such products 
have a singularity. In this paper we are considering interacting fields so it is not known if a singularity would 
arise for such products of operators evaluated at the same point. Heisenberg in his investigations of a quantized 
nonlinear spinor field repeatedly underscored that in a quantum field theory with strong interaction the singu- 
larities of propagators can be essentially smoothed out. In this situation the problem of the renormalization of 
nonrcnormalizable QFT will be crucially changed. This problem appears when we have the product of propa- 
gators of some quanta, but the nonrenormalizable theories have to be considered on the nonperturbative level 
which probably gives rise to smoothed out, nonsingular Green's functions. Thus one might suppose that the 
nonperturbative quantization methods may resolve the problem of the quantization of nonrenomalizablc QFT. 
Physically it is hypothesized that there are situations in interacting field theories where these singularities do 
not occur (e.g. for flux tubes in Abelian or non-Abelian theory quantities like the "electric" field inside the 
tube, (£%) < oo, and energy density e(x) — ((£~) 2 ) < oo are nonsingular). Here we take as an assumption that 
such singularities do not occur. 



3 The calculation of 2 and 4-points Green's functions 

The quantization of equations (|2.6|) - (|2.7|) in the spirit of section [5] evidently is very hard and deriving exact 
results is probably impossible. In order to do some calculations we give an approximate method which leads to 
the 2 and 4-points Green's functions only. In order to derive the equations describing the quantized field we 
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(3.1) 



average the Lagrangian over a quantum state | Q) 

Q £sU(3)(x) = (£sU(3)) = 

\gf BCD ({d^{x) - d v A%(x)) A c »(x)A D »(x)) + 

^ g 2 f BC lDl f BC 2 D 2 (A^{X)A^{X)A C ^{X)A D ^{X)) . 

Now we will detail the kind of physical situations we wish to describe. The model given here is similar to 
stationary turbulence when there are time dependent fluctuations in any point of the liquid but all averaged 
quantities are time independent. For a QFT this means that all Green's functions are time independent and 
there is a correlation between quantum fields in different points at one moment 

{A*H&)-A%ixZ)) + 0, (3.2) 
h = ■ ■ ■ = t n , (3.3) 
r! ? ■■■^r n . (3.4) 

In linear and perturbative QFT this is not the case because the interaction is carried by quanta which move 
with a speed less than or equal to the speed of light. In these theories the correlation between quantum fields 
in different points at the same time is zero. 

= o, (3.5) 

h = ••• = *„, (3.6) 
ri ? ■■■ + r n . (3.7) 

In this sense one can say that nonperturbative QFT in some physical situations is very close to turbulence, i.e. 
in nonperturbative QFT there may exist extended objects where quantized fields at all points are correlated 
between themselves (example from QCD are flux tubes and glueballs). Such objects fall into two categories 

1 . The averaged value of all quantized fields are zero 

(A*(x))=0. (3.8) 

But the square of these fields are nonzero 

\{A*{x)Y)^Q (3.9) 

2. Some components of the quantized fields are nonzero and some zero, and the square of some components 
is nonzero 

(-4^ 0*0) + forsome ^€1,2, ■■-,8 (3.10) 
(A® 2 (x)) = for remaining B 2 e 1, 2, ■ ■ ■ , 8 (3.11) 

(Af 3 {x)) 2 } ^ forsome S 3 el,2,---,8 (3.12) 

The most natural case is when A^ 1 belongs to a small subgroup of SU(3) gauge group (for example, to 
SU(2)) and A^ 2 are the coset components and B 3 = B 2 . For example, B\ = 1, 2, 3 and B 2 = 4, 5, 6, 7, 8. 

In the first case the quantized fields are in a completely disordered phase. In the second case one has both 
ordered and disordered phases. In this paper we present some calculations which give illustrate each case. The 
first case of the completely disordered phase will be illustrated by a glueball configuration; the second case of 
both ordered and disordered phases will be illustrated by a flux tube solution. 



4 The first case: completely disordered phase — glueball 

4.1 The vector to scalar field replacement approximation - — > (fi B 

In any quantum field theory the Green's functions contain the full information about the quantized fields. In 
this section we will present equations which describe only the 2 and 4-points Green's functions for SU(3) gauge 
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theory (i.e. QCD) using the approximation where all the components of the gauge potential are in a disordered 
phase. To begin we average the SU(3) Lagrangian using some approximate expressions for the 2 and 4-points 
Green's functions. The SU(3) Lagrangian is 



C 



SU(3) 



- 



(4.1) 



where T B V = d^A 8 - d v A B + gf BCD A^A^ is the field strength operator; B, C, D = 1, . . . , 8 are the SU(3) 
color indices; g is the coupling constant; f BCD are the structure constants for the SU(3) gauge group; A B is the 
gauge potential operator. In order to derive equations describing the quantized field we average the Lagrangian 
over a quantum state | Q) 



Q 



c 



Q 



c 



8 A 



d^Ai 



d v A B » 



\gf BCD {(d»A B -d v A B )A c »A D » 

1 2 f BC lDl f BC 2 D 2 (jCifih 

This expression contains the following 2, 3 and 4-points Green's functions: 
Defining the 2-point Green's function 

5^(x,y), 



(4.2) 



'A B (x)A^(y) 
the first term on the r.h.s of equation (|4.2|l becomes 



t A B {x)) (d»A B »{x)) 



= d x ^ (a b (x)) (A B »{y) 



where we denote as x and y^ as y. For simplicity we consider the case with x° = y° 
function we will make the "one-function approximation" 

rABr \ j-ACD sBCE iD i \iEi \ 

Gap (X,y) ~ -Va/3j f <P (X)4> (y) 



(dAf), {(dA)A 2 ) and {(A) 4 

(4.3) 

(4.4) 

For this Green's 



(4.5) 



where <j> A [x) is the scalar field which describes the 2-point Green's function. Other authors have used a similar 
technique of replacing a vector, gauge field, by some combination of scalar fields. For example Corrigan and 
Fairlie [IHj and Wilczek [HZj were able to transform the Yang-Mills equations into those of a massless scalar A^> 4 
theory by writing the SU(2) gauge fields in terms of of some combination of effective scalar fields. Physically 
this approximation means that the quantum properties of the field A B can be approximately described by a 
scalar field (j> B (x), i.e. in this approximation the Lorentz character and Lorentz index /i of the vector gauge 
fields are treated in a simple way through H4.5JI . Taking into account this approximation we have 



(d,A B ) (d»A Bv )) = ~ifJ BAC f BAD (d^ c ) (d^ D ) = -12 (d^ A ) (d^ A ) 



and 



[d,A B ) (d v A B »)) = -3 (d^ A ) (d^ A ) , 
Next we make the assumption that the odd Green's functions average to zero 

(A B (x)A%(y)A D ;(z)) = 0. 

This gives 

((d»A B (xj) A c p (x)A^{x)) = d x , (A B (x)A c (y)A^(z) 

which finally leads to 

~ ' 







d^A B - d v A B ) A c »A Dv 
Last for the quartic term on the r.h.s of equation (|4.1|) we make the approximation that 

[A B {x)A%{y)A°{z)Af{u)) « (A B (x)A°(y)) (A°(z)Af(u)) + 

'A B (x)A^(z)) (A%(y)A?(u)) + (a b (x)A*(u)) (A c (y)A°(z) 



(4.6) 
(4.7) 

(4.8) 
(4.9) 
(4.10) 

(4.11) 



Nonpcrturbative calculational method in quantum held theory 



7 



This assumption is that the 4-point Greens function is simply the sum of products of the 2-point Green's 
function. In this approximation the l.h.s of l|4.11f) is 



A B ^x)A^{x)A D ^x)A Rv { X y i 

Al, 2 ;(P 12 ,Q 12 )(/ B£lP V C£lQ ^^^^ (412) 
Al, 2 ;(P, 2 ,Q, 2 ) (f***f°**4* (X)4fl* (*)) (f CE ^f RE ^4> P >{x)^{ X )) + 
Al, 2 ;(P 12 ,Q 12 )(/ BBlPl / ii£lQl Pl W^ 1 (-)) (f^f^^^ix))^ 



We have taken the spacetime variables of the scalar field, </>, such that x = y = z = u. In the above expression 
^i,2;(Pi 2 ,Qi 2) are some parameters depending on the values of the indices Pi, 2, Qi,2 

!Ai, if all indices Pi, 2, Qi,2 = 1, 2, 3, 
A 2 , if all indices Pi, 2 , Qi, 2 = 4, 5,6, 7, 8, (4.13) 
1, otherwise 

where Ai, 2 are constants. Equation (|4.13|) essentially introduces three different coupling strengths (Ai, A 2 , and 
1) in place of the original single coupling. This step is necessary in order to insure that the solutions we find 
are at least quasi-stable. This will be explained in more detail when we discuss the numerical solutions of this 
system in the next section. From the definition of the index Ai^rpj 2 ,q 1 2 ) one can see that this approximate 
quantization procedure acts differently for the scalar field components which belong to the small subgroup 
517(2) S SU(S) versus the coset SU(3)/SU(2). The first term in eq. l|4"l"2"|) can be written as 

tfBExPi fCEiQi^Pi^Qi^ ^DE 2 P 2 jRE 2 Q 2 ^P 3 ^Q 2 ) 

A, (/ S£lQ / CBl V a fc ) (f DE2C f RE2d (t) c (t> d ) + (4-14) 
A 2 (/ BE i m / c ^^ m f) (f DB *P f jfl) + (other terms) 

A straightforward calculation using the antisymmetry property of the structure constants shows that 

j ABC jyADR (jBPiPi jCBiQi^Pi ^Qi\ ^DE 2 P 2 jRE 2 Q 2 ^2^)2^ _ 0) (4.15) 
f ABC f ADR^BE 1 a f CE 1 b r(j) b^ f DE 2 c f RE 2 d (j) c^ = ^ (41g) 
jABC fADR ^fBEnn fCEm^m^n^jDE-ip fRE 2q( jjp^q^ _ q (4-17) 



which consequently yields 

^MP lM f AB °f ADR {fBE^ f CE lQ ^P^ Ql) 
( f DE 2 P 2f RE 2 Q 2< pP 2< p Q2)ll ^ =Q 

A similar calculation for the second term in eq. I|4.12[l shows that 



(4.18) 



Consequently 



j ABC f ADR ( f BE lPl f^Q^P.^Q,^ (4 1Q) 
( f CE 2 P 2f RE 2 Q 2(t) P 2{x) ^ 2{x)) = W {r(j) a + r())m) ^ 

f ABC f ADR^B El a f DE lbr(j) b^ f CE 2 c f RE 2 d^^ = £7 ( ^ a)2 ^ ^ ^ 

97 

j ABC j ADR. ^jBEim jDEin^m ^n^j (^jC E 2 p jRE 2 q ^ = f± (^™0™) 2 , (4.21) 



K2;(P lM f ABC f ADR (f^f^Q^ P ^ Ql ) (f^P 2f RE 2 Q 2 ^, 



27 o 27 o 27 (4.22) 

-Ai {rrf + (0"v m ) + -r (<w (<T<T) ■ 

o 4 



Analogously the third term in eq. (|4.12() can be shown to give 



\2;(?iMf ABC f ADR (f BE ^f RE ^^^) ( f CE 2 P 2f DE 2 Q 2(/) P 2(j)Q2) = 

27 2 27 2 27 ( 4 - 23 ) 

T 



-Ai (4> a 4> a f + — -A 2 >^ m ^ 2 ' " <•■"•-"> 
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Combining these results gives the following expression for the quartic term 

/ \ 81 81 

jArbjAcd /aWA c »A d v) = (<j> a <t> a f + ttA 2 (0 m m ) 2 + 81 (0 a a ) (<j> m <p m ) . (4.24) 



Bringing together all the results gives the following effective Lagrangian which describes the 2 and 4-points 
Green's functions 



yAr (rrf + yA 2 {<T<Tf + 81 W) (<F<T) 



(4.25) 



This Lagrangian can be put in the standard form by making the following redefinitions: <f> a — > 2(f> a /(3g) and 
Ai,2 — > Ai j2 /2 which yields 

y£ e // - (9^ A ) + y (<T<T) 2 + y (0 m m ) 2 + (W (0™0 m ) . (4.26) 

By splitting up the kinetic energy term into parts for cf> a and cf> m one sees that this Lagrangian is of the form of 
two, massless scalar fields with A0 4 self interaction and an interaction term between the two scalar field (i.e. the 
last term in iET^o) ). Now in [SHI [t 

was shown that radiative corrections could introduce a symmetry breaking 
(i.e. negative) mass term into a scalar Lagrangian similar to that in (|4.2(il) . This effect is called dimensional 
transmutation. For the pure scalar case one can not rigourously justify a radiatively generated symmetry 
breaking term since the scale at which the symmetry breaking occurs lies outside the region where pertubation 
theory is valid. Nevertheless it was postulated that a nonpertbative calculation would yield a similar negative 
mass, symmetry breaking term. In the present paper we will also make this assumption and add negative mass 
terms to the two scalar fields so that the Lagrangian can be written as 

4 2 4 4 ( 4 27 ) 

-f rv" 1 - ffffl 2 - -j (co 2 + W) (<r<n 

where (f> A = (<pQ, 4>™) are some constants. The mass terms come from the cross terms of the second and fourth 
terms above. Next we drop the constant terms in the Lagrangian which gives 

^£ e// = -\ {d^ A f + T irr - <t> a M? + y [<t> m <t> m - effl 2 + (rr) (<r<n . (4.28) 

In the next subsection we will investigate numerically if this Lagrangian leads to a solution which is well behaved 
in some finite range, and which therefore could be viewed as a glueball. As a first step we write down the field 
equations that give an effective scalar field description of QCD 

= [2ff + Ai (rr - «)] , (4.29) 

d^f" = -f[2ff + A 2 (ff-ffff)]. (4.30) 

In ref. |59| a similar coupled scalar field system was investigated and shown to have soliton solutions for some 
specific choice for the potential terms. In the next subsection it will be shown that the Lagrangian in 14.28fl has 
solutions are well behaved in some finite range. We will interpret these solutions as glueballs since the original 
starting Lagrangian was the pure SU(3) gauge theory. The importance of the quantization procedure that led 
to the Lagrangian in (|4.28|) is evident. For classical, pure SU(3) theory one has no finite energy solutions - there 
are no classical glueballs. Our approximate, nonperturbative quantization method transformed the original pure 
Yang-Mills into an effective interacting scalar Lagrangian of (|4.28|) . In the next section solutions which have 
a soliton-like character in some finite range are investigated, and it is suggested that this solutions may give a 
model for glueballs. 

4.2 Soliton/Glueball solution for the (f) A field 

We begin our investigation of the field equations (|4.29|1 l|4.30(l by taking the following ansatz for the scalar fields 

^(r) = a =1,2,3, (4.31) 

0»( r ) = ZM 5 m = 4,5,6,7,8. (4.32) 
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This ansatz implies that the components a and the components m have different behavior from one another. 
Substituting equations (|4.31[1 (|4.32[) into equations l|4.29[l (|4.30[) gives the following coupled, nonlinear equations 

0" + V = [f 2 + Ai (0 2 — m 2 )] , (4.33) 

/" + ;/' = /[0 2 +A 2 (/ 2 - M 2 )] (4-34) 

The primes denote differentiation with respect to r, we have redefined Ai, 2 /2 — » Ai, 2 , and 20g</>o = ™ 2 anc ^ 
20™ 0™ — /j 2 ; m,fi are constants which will be calculated by solving equations (|4.33|l and (|4.34(l . We have 
not found analytical solutions to these equations and so we investigated them numerically. The preliminary 
numerical investigation showed that these equations do not have regular solutions, f(x),(f>(x), for an arbitrary 
choice of the m, jj, parameters. However, certain discrete values of these parameters appear to give solutions 
which are well behaved at least for some finite range. Thus we will solve equations l|4.33[) l|4.34|l as a nonlinear 
eigenvalue problem for eigenstates (f>(x),f(x) and eigenvalues m,fi, i.e. we will calculate the parameters m,/x 
such that the solutions 0(r) and f(r) are regular functions for some finite range. 

The solutions of l|4.34|l will depend on the initial conditions of the ansatz functions (0(0), /(0), 0'(O), /'(0)) 
and on A1.2. We will rescale the solutions by dividing equations l|4.33l) (|4.34|l by 3 (O), and then introducing the 
dimensionless radius x — r0(O). Finally by redefining (f>(x) / 0(0) — > <j>(x), f(x)/(f>(0) — ► f(x) and m/<f)(0) — ► m, 
/j,/ <fi(0) — > /1 we arrive at the rescaled equations 

0" + V = 0[/ 2 + A 1 (0 2 -m 2 )], (4.35) 

/"+-/' = /[0 2 + A 2 (/ 2 - M 2 )]. (4.36) 

where now the primes denote differentiation with respect to x. We will look for solutions to l|4.35|l 1|4.36|) which 
are finite at the origin, and which approach finite values as x — » 00. In particular we would like the asymptotic 
values to approach (f>(oo) = m and f(oo) = 0. We have not been able to numerically integrate the equations 
out to arbitrary x, but we have found solutions with 4>(xf) — Tn and J(xf) = 0, where xf is some finite x 
coordinate. To implement the numerical solution of (|4.35|) (|4.3()|) we used the following boundary conditions 

0(0) = 1, <£'(0)=0, (4.37) 
/(0) = fo, /'(0) = 0. (4.38) 

These boundary conditions will give solutions which have soliton-like profiles for the region < x < xf- We 
now rewrite equation (|4.36(l in the following form 

- (/" + V) + fVeff = (A 2M 2 ) / (4.39) 
where we have introduced an effective potential 

Veff - {<t> 2 + A2/ 2 ) • (4.40) 

Thus equation l|4.39|l is now in the form of the Schrodinger equation. It will have regular solutions provided 
V e ff has a deep enough "potential well" . Taking into account the desired boundary conditions of 4>(x) and f(x) 
for x — > xf (i-e. (f>(xp) — > m and f(xp) — > 0) we want V e ff(xp) > V e ff(0) or m 2 > 1 + A2/0 • We would like to 
set xf = 00, but numerically we have only been able to integrate the equations l|4.35Jl (|4.36|l out to some finite 
xf- If m 2 is too close to 1 + A2/0 the "potential well" may be too shallow and no solution exist. However, 
assuming m 2 is large enough a regular solution will exist and the "energy levels", A 2 /i 2 , will be quantized. 

4.3 Numerical solution 

We will use the following numerical method for solving equations l|4.35|l l|4.36() : input some arbitrary zeroth 
order approximation for the function f(x) (which is denoted fo(x)) and solve equation l|4.35(l for the zeroth 
order (f> (</>o) 

K + -Jo - 0o [/o 2 + Ai (4>l - m 2 )} (4.41) 

where mo is the zeroth order approximation for the parameter m. The boundary conditions are given in l|4.37|l . 
This regular solution (f>o(x) is then used to calculate the next order f(x) ansatz function. Taking go(x) and 
substituting it into equation l|4.36|l we now have an equation for the next order in f(x). 

/r + ^ = A[0^ + A 2 (/ 2 -M 2 )] (4.42) 
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with the boundary conditions of equation l|4.38[) (for the present numerical calculations we take /o(0) = V0.6). 
Numerically solving l|4.42J) then gives the first approximation fi(x) which is then substituted into equation 



(4.43) 



This equation is then solved for the first approximation in <pi (x) which is then put back into the next approximate 
equation for f(x) and so on i.e. 



and 



< + ~J\ = & [fli + Ai {<t>1 ~ ml)] 



x 



(4.44) 
(4.45) 



At each step we find values mf and /sf which approach the true eigenvalues values m* and n* . 



4.3.1 A more detailed description of the numerical calculations 

For the numerical solution of equation 14.41fl we make an initial arbitrary choice of the zeroth order approxi- 
mation for f(x) as 

cosh § 



The reason for this choice is that the cosh function in (|4.4(i|) has a soliton-like shape, and we are looking for 
soliton-like solutions. Typical solutions for arbitrary values of to,q are presented in Fig^ From the figure one 



=1.837435 




Figure 1: The typical singular solution for equa- 
tion (OTJ}. Ai = 0.1. 



H!=1- 492103 



10 15 



Figure 2: The typical singular solution for equa- 
tion g32l- A 2 = 1.0. 



sees that when mo < TOg (toq is the unknown parameter which is conjectured to give a regular solution) the 
solution <j)Q (x) is singular. Near to the singularity the equation Ij4.41|l has the form 



and the following approximate solution 



0o(x) 



1 



Ai xq — x 



(4.47) 



(4.48) 



where xq is some constant depending on too. On the other hand for too > Toq a typical solution is shown in 
FigGl The corresponding asymptotical equation is 



2 

^(x) + ~<t>' a w - (Aito 2 ) 0o 



(4.49) 



which has the following approximate solution 



4> (x) s» <j) c 



sin ( x \J\\m^ + a ) 



(4.50) 
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where <j)oo and a are constants. We conjecture that there may exist some value, mj, which yields an exceptional, 
non-singular solution. Within some numerical accuracy such a solution is presented in Fig0] For this value 
the equation (|4.41|) has the following asymptotical behavior 



+ 2Ai (m*) 2 (0o - ml) 



(4.51) 



and the corresponding asymptotical solution is 



4> Q (x) ml + (3- 



x\J 2Ai (r. 



(4.52) 



where (ft^ and [3 are some constants. 

The next step is finding the first approximation for the fi (x) function. The equation is 



(4.53) 



From the previous calculations we use the exceptional regular solution 4>o(x) with the asymptotical behavior 
(I4.52|) . Then the numerical investigation shows that for an arbitrary /i there are two different singular solutions 
for fi (x) which are presented in Fig|21 As with equation (|4.47|) the singular behavior of the function fi (x) is 





Figure 3: The iterative functions /i, 2, 3, 4(2;)- 



Figure 4: The iterative functions 4> 1,2,3, a{x). 



given by 



A(*) 
hi?) 



1 



A 2 x - x 
2~ 1 



A 2 x - x 



by (Mi < /i*, 
by fix > n\. 



This again may indicate that there is some intermediate special value for [i\ (i.e. [i\ 
solution, f*(x), is regular and has the following asymptotical behavior 



/*(*)«/«> 



(4.54) 
(4.55) 

fi\) for which the 
(4.56) 



where /qq is some parameter. The next step is substituting the first approximation f*(x) into equation 14.3511 to 
find the regular, exceptional solution 4>\(x) with eigenvalue m\\ then 4>\(x) is substituted into equation l|4.36(l 
to find the regular exceptional solution /^(x) with eigenvalue fi = fii, and so on. 

The result of these calculations is presented in Fig's. 00] and Tabled We see that there is the convergence 
4>*{x) — > <P*(x), f*(x) — * f*(x), m* — > m* and [i* — > f*(x),(j)*(x) are the eigenstates and m*,fi* are 
eigenvalues of the nonlinear eigenvalue problem (|4.33H 14.340 . 

From the figures it can be seen that the numerical integration of the coupled, nonlinear equations was only 
carried out to a finite x (roughly x — 20). Beyond this it became increasingly difficult to adjust the mass 
values, m* , fi* in order to obtain the desired soliton-like behavior. At this stage it is not certain whether the 
solution can be made to be well behaved out to arbitrary x. However, the analogy between our system and 
the Schrodinger Equation given in l|4.39|l indicates that there should be some parameters which give regular 
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i 


i 


2 


3 


4 




1.590537877.. . 


1.579111774.. . 


1.565224661... 


1.5679008386... 




1.222955... 


1.2287.. . 


1.22902.. . 


1.229019.. . 



Table 1: The iterative parameters m* and /i*. 



solutions out to x — oo. However based on the numerical results it is at present not certain as to the asymptotic 
behavior of the ansatz functions. 

This numerical method was tested by applying it to the standard, analytic soliton solution of A</> 4 theory in 
1 spatial dimensions. Details of this can be found in Appendix FBI 

Finally, we remark that according to Derrick's Theorem |61| our solution, if it turns out to be regular for 
< x < oo, it will not be stable. For pure scalar field theories as in (|4.28|) absolutely stable solutions exist only 
in 1 or 2 spatial dimensions. Thus our solution is at best quasi-stable. This is in fact the reason for choosing 
different couplings (i.e. Ai,A2, 1) in (|4.13|1 . This resulted in the potential for the interacting scalar fields of 
(|4.28l) having global and local minima. Our solution is in one of the local minima. Thus although it is not 
absolutely stable it should be quasi-stable so long as there is a large enough barrier between the local and global 
minima. The can be accomplished by requiring that Ai and A2 be significantly different. This however raises 
the problem of why one should develop two different scales from a theory that originally only had a single scale. 
We will leave this open question for future work. 

The final remark is that we can add the constant A2 (t/'cT'/'o 1 ) 2 (which have not any influence on the field 
equations l|4.29|l (|4.30|l ') to the Lagrangian l|4.28|l we will have the finite energy of glueball. The profiles of the 
enrgy with respect to A1.2 and with different (fio are presented in fig's 15161 



-■-^ = 1.0 -1-^=1.0 

-•-^ = 1.5 -t-^aLS 




Figure 5: The profile of the glueball energy with Figure 6: The profile of the glueball energy with 

respect to Ai. respect to A 2 . 



5 The second case: ordered and disordered phases — flux tube. 
5.1 Reduction to a small subgroup 

In this section we will describe the second case where one has both an ordered and disordered phase. We 
reiterate that we have not found an exact analytical method for solving the full equations (|2.5(l (|2.6|l for all 
the Green's functions. Our basic approach in this case is to give some physically reasonable scheme for cutting 
off the infinite set of equations for the Green's functions. In addition our approximate calculations will involve 
the decomposition of the initial gauge group into a smaller gauge group: SU(3) — ► SU(2). Physically we will 
find that this reduction splits the initial degrees of freedom (the SU(3) gauge potential components) into SU(2) 
and coset components. The SU(2) components will represent the ordered phase while the coset components, 
which again will be represented via an effective scalar field, will represent the disordered phase. After these 
approximations one can perform analytical calculations which suggest the formation of a color electric flux 
tube. This is of interest since an important feature of the confinement phenomenon in the dual superconducting 
picture of the QCD vacuum is the formation of electric flux tubes between the confined quarks. 

First we define the reduction of SU(3) to SU(2) following the conventions of ref. [21] • Starting with the 
SU(3) gauge group with generators T B , we define the SU(3) gauge fields, = A B T B . SU(2) is a subgroup of 
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SU(3) and SU(3)/SU(2) is a coset. Then the gauge field A^ can be decomposed as 

A? = A B T B = a«T a + A™T m , (5.1) 
al e SU(2) and A™ G SU(3)/SU(2) (5.2) 

where the indices a,b,c. .. belongs to the subgroup SU(2) and m, n, . . .to the coset SU(3)/SU(2); B are SU(3) 
indices. Based on this the field strength can be decomposed as 



T* V T B =T%T a +3^ v T m (5.3) 

where 



^ = h%, + *% eSU(2), (5.4) 

K» = d^a a u -d v al + ge abc alat e SU(2), (5.5) 

*%, = 9f amn A™A n v &SU{2), (5.6) 

^Tu = F™+G™ £SU(3)/SU{2), (5.7) 

F™ = d lt A?-d v A'£+gr mp A»Al e SU (3) / SU (2) , (5.8) 

= 9 r nb ( An A- A >\) £SU(3)/SU(2) (5.9) 

where f ABC are the structure constants of SU(3), e abc — f abc are the structure constants of SU(2). The SU(3) 
Yang-Mills field equations can be decomposed as 

d v {h a ^ + ^ v ) = -gf amn A^{F n ' iV + G n i u '), (5.10) 
D v (F m ^ + G m H = -gf mnb [A n v (h b ^ + - a h v (F n ^ + G n ^)] (5.11) 

where d v [- ■ - ] a — d v [- ■ - ] a + f abc a b [- ■ - ] c is the covariant derivative on the subgroup SU(2) and D v [- ■ -] m = 
d u [- ■ ■ ] m + f mnp A™[- ■ ■ ] p . The Heisenberg quantization procedure gives us 

dv {h a ^ + S "") = -gf amn A™ (F ntlv + G nA1I/ ) , (5.12) 
D v (P m ^ + G m ^) = -gf mnb [JSJJ (h bt * u + - a b u {P n ^ + G n ^)] (5.13) 

5.2 Basic assumptions 

It is evident that an exact quantization is impossible for eqs. (|5.12|) Q5.13[) . Thus we look for some simplification 
in order to obtain equations which can be analyzed. Our basic aim is to cut off the set of infinite equations using 
some simplifying assumptions. For this purpose we will propose simple but physically reasonable ansatzen for the 

2 and 4-points Green's functions - (A^(y)A^(x)) 1 /a%(x)a b 3 (y)A^(z)A^(u)^ and (A%(x)A%(y)A*(z)Al(u)\ - 

which are involved in averaged Lagrangian l|4.2|) . As was mentioned earlier we assume that there are two phases 

1. The gauge field components (a=l,2,3 £ SU(2)) belonging to the small subgroup SU(2) are in an 
ordered phase. Mathematically this means that 

= K{x)) cl . (5.14) 

The subscript means that this is a classical field. Thus we are treating these components as effectively 
classical gauge fields in the first approximation. 

2. The gauge field components A™ (m=4,5, ... , 8 and A™ <E SU{3)/SU(2)) belonging to the coset 
SU(3)/SU(2) are in a disordered phase (i.e. they form a condensate), but have non-zero energy. In 
mathematical terms this means that 

(A™(x)) = 0, but (A™(x)A n v (x)) + 0. (5.15) 

Later we will postulate a specific, and physically reasonable form for the non-zero term. 

3. There is no correlation between the ordered (classical) and disordered (quantum) phases 

(/KMO) = /«) (g(A™)) (5.16) 

Later we will give a specific form for this for the "mixed" 4-point Green's function, 
(a a a (x)a b (y)A^(z)A^(u)) 
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5.3 Derivation of an effective Lagrangian 

Our quantization procedure will deviate from the Heisenberg method in that we will take the expectation of 
the Lagrangian rather than for the equations of motions. Thus we will obtain an effective Lagrangian rather 
than approximate equations of motion. The Lagrangian we obtain from the original SU(3) pure gauge theory 
will turn out to be an effective SU(2) Yang-Mills-Higgs system which has monopole solutions. The averaged 
Lagrangian is 

c = -\ {t%t a ^) = -~ «^vr^) + (r™r m n) ( 5 - 17 ) 

here and T m » v are defined by equations l(0 )l -lfO |l . 

5.3.1 Calculation of (T^T^) 

We begin by calculating the first term on the r.h.s. of equation (|5.17() 

{T%T a ^) = (h« u h a ^) + 2 (ht a ,$*v) + ■ (5.18) 

Immediately we see that the first term on the r.h.s. of this equation is the SU(2) Lagrangian as we assume that 
a° and /i" are effectively classical quantities and consequently 

(h^h^) « h%h a ^. (5.19) 

The second term in equation H5.18H is 

/ jl a rf.aLLis\ ramn //a a a a\ Amu Anu\ , rabc ramn /be Amu, atiu\ (vr on\ 

{h^u® )=9J \{Op,a u - c^aj A ) + gf f {a^a^A ^ A ). (5.20) 
Using assumptions 1 and 3 from the previous section these terms become 

(a a a (x)A™(y)A»(z)) = a a a (x) {A™{y)A n v {z)) = n^a a a (x)G mn (y,z) (5.21) 

and 

(a a a {x)a b p {yW;{z)Al{u)) = a a a {x)a b {y) Vllv g mn {z, u) (5.22) 

The function Q mn (x,y) is the 2-point correlator (Green's function) for the disordered phase. Because of the 
bosonic character of the coset gauge fields G mn (x,y) must be symmetric under exchange of these fields. Also 
by assumption 2 of the last section this expectation should be non-zero. We take the form for this 2-point 
correlator to be 

(A™(y)A?(x)) = ~ l -^r pb r pc 4>\vW^) = -V^Q mn (y, x) (5.23) 

with 

mn (y, x) = \r pb r pc <t>\yW{x) (5.24) 

here <f> a is a real SU(2) triplet scalar fields. Thus we have replaced the coset gauge fields by an effective scalar 
field, which will be the scalar field in our effective SU(2)-scalar system. Aside from the factor of — ^ this is 
similar to the approximation made for the 2-point correlator in the last section in l|4.5|) . The factor of — ^ is 
introduced so that the effective scalar field, tj>, will have the correct coefficient for the kinetic energy term. If we 
were to take the scalar field to be constant (<fi a (x) « const.) then (|5.23l) and (|5.24|) would represent an effective 
mass-like, condensation term of the coset gauge fields. With this we find that the middle term vanishes 



(h%^ v ) = grf" ( f amn ~ S mn (x, x)+ 

(5.25) 

gf amn f abc a b Jx)at(x)g mn (x,x)) =0 



The last term which is quartic in the coset gauge fields will be considered at the end. Up to this point the 
SU(2) part of the Lagrangian is 

(F^J 70 ^) = (h a ^ v h a ^) + g 2 f an Pf an 'p' (A$A*A n 'i>A ffu } . (5.26) 
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Analogously 

2 

Next we consider terms with three coset fields. The term involving the derivative is 

f mnp ((d y »A™ (y))A^ {x)A pv (z)) (5.36) 

Since the gauge fields must be symmetric under exchange, and because of the antisymmetry of the of f mn P this 
term vanishes. Next the terms involving three coset fields and one SU{2) field we will approximate as 



jbrnn' j^mnp^b ^ A n ' j^nf 1 j\P vs ^ ~ — jbrnn' jirnnp ^b ( I^J^'^jJ^'^A? 1 '} 



(5.37) 



(A™ A nii ){A pv ) + {A£ A p,y )(A nt *) 
By the second assumption in the previous section, (AV? (x)) — 0, this term also vanishes. Thus 



2 

(T™T mllv ) = -2d ti (j) a d ll ct> a + 2ge abc a a ^d t ,cj) b (t) c - ?-e abc e ab ' c a b (t> c a b • ' 



-fLaffiaWj," + gZpnnpjmn'p' /^^/d/^ = (5.38) 

-2 (d^ 11 + |e abc a^ c ) - yo^'a^f + g 2 f mn Pf mn 'p' ^A p A n '^A p ' 1 



The full averaged Lagrangian is 

-\ (^ A n = - \k^ v + \ + f ^ &c «: 

- g 2 f An Pf An 'p' (A^APA n '^A p '^ . 
where we have collected the quartic terms from eq. (J5.26JI and (|5.38(l together into 

jAnp j-An'p' / A n j\P A n ' H J[p' " 



(5.39) 



5.3.3 The quartic term 

In this section we show that the quartic term - f An P f An v I^A , ^A v v A n A p - from cq. Ij5.51|l becomes an 

effective A^ 4 interaction term for the effective scalar field introduced in eq. (|5.24(l . Just as in eq. I|5.24|l where 
a quadratic gauge field expression was replaced by a quadratic effective scalar field expression, here we replace 
the quartic gauge field term by a quartic scalar field term 

(A^(x)A^y)A^z)Al(u)) = (e™Z^ v + K^d^ 0v + E^t,^) 

here E" l ^ d , E^ p ^ d , E nl ^ p d are constants. Because of the bosonic character of the gauge fields in H5.40J) the 
indices of these constants in conjunction with the indices of the r) a ps must reflect symmetry under exchange of 
the fields. The simplest choice that satisfies this requirement is 



{A™{x)A%y)Al(z)Al{u)) = (S mn S^r, a ^ u + S mp S n ^ V0u + <T^ M /3m) 

e abcd r(x)<t> b (y)<f> c (z)<l> d (u) 



(5-41) 



This choice of taking the constants from eq. (|5.40|) to be products of Kronecker deltas and fixing a = b = c = d 
for the lower indices, satisfies the bosonic character requirement for the gauge fields, and is equivalent to the 
reduction used for the quartic term in (SHIED!- Evaluating eq. (|5.41(l at one spacetime point {i.e. x = y = z = u) 
and contracting the indices to conform to quartic term in eq. I|5.51|l gives 

(AZ(x)Aftx)A*{x)Al{x)) = (*™a»VTr + S mp S n ^X + 5 m % np rf^) 

9 (5.42) 
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where the constants e a bcd are chosen that at the point x = y = z = u 

eabcd = ~ {SabScd + SacSbd + SadSbc) ■ (5.43) 

This expression can be further simplified to 

(A^A n v A^A qv ) = (45 mn 6 pq + 16S mp S nq + 45 mg 5 np ) (4> a (x)4> a (x)) 2 . (5.44) 
Substituting this into the original quartic term of eq. I|5.51|l yields 

jAnpjAn'p' ^j^n jg> j^n' 'n j^t 'v\^ _ 
^jAnnjAn'n' + ^jAnpfAnp + ^AnpjApr^ ( x )f = (5.45) 

12f Anp f Anv (4> a (x)4> a (x)f 

where the antisymmetry property of the structure constants has been used. Using the explicit expression for 
the structure constants (/ 123 = 1; / 147 = / 246 = / 257 = / 345 - / 516 = / 637 = |; / 458 = ,/ 678 = & plus those 
related to these by permutations), and recalling that the index A runs from 1 — 8 while the indices n,p run 
from 4 — 8 one can show that f An Pf An P = 12 (/ 123 and related constants do not contribute to this expression). 
Combining these results transforms the quartic term in eq. (j5.51(l as 

l g 2 f An Pf An' P ' ( A « A P A n'* A P'»} = 36<? 2 (x)4>*{x)f = j (x) 0« {x)f (5.46) 

One can show that the simplification (|4.11|1 is also applicable here and will lead to the same result. This has 
transformed the quartic gauge field term of the coset fields into a quartic interaction term for the effective scalar 
field. Substituting this result back into the averaged Lagrangian of eq. (|5.51(l we find 

-- (^ A n = -\h^h^ + \ (d,<p a + f e abc a b ^ c ) 2 + ^ a ^W- 

4 4 I \ I \ ( 5 ' 47 ) 

-{r{xw{x)f 

The original pure SU(3) gauge theory has been transformed into an SU(2) gauge theory coupled to an effective 
triplet scalar field. This is similar to the Georgi-Glashow Lagrangian except for the presence of the term 
^-a b 1 (j> b a c,1 (j> c and the absence of a negative mass term for <f> a of the form m 2 (j) a (x)(j) a (x). 

The Georgi-Glashow Lagrangian is known to have topological monopole solutions |63j which have the form 

r = ^-, < = 0; at= ealbXb[1 - h(r)] (5.48) 

where /(r) and h(r) are functions determined by the field equations. For this form of the scalar and SU(2) 
gauge fields the term, ^-a* <fr b a cfJ "(f) c , vanishes from the Lagrangian in eq. i|5.47[) by the antisymmetry of e a ib and 
the symmetry of x a x b . Thus for the monopole ansatz of eq. (|5.48|l the Lagrangian in (|5.47|) becomes equivalent 
to the Georgi-Glashow Lagrangian minus only the mass term for the scalar field. 

In the present work we simply postulate that the effective scalar field develops a negative mass term of the 
form — ^(4> a <fi a ) which is added by hand to the Lagrangian of Ij5.47|l to yield 

2 (5.49) 

The scalar field now has the standard symmetry breaking form and this effective Lagrangian has finite energy 
't Hooft-Polyakov solutions (the last term should not alter the monopole construction since it vanishes under 
the ansatz in l|5.48|l 'l. 

At the end of this section we have to note the following problem. The covariant deravitives D^af, — 
d^a^ + ge abc a b ^al and Z? M <^ a = d^(j) a + f e ahc a^0 c are different. One can avoid this problem if we suppose that 
there is a correlation between ordered phase and disordered phase A™ 

«i • • ■ <*£ (x n )A™(y)A}(z)) = k n a% {x x ) ■ ■ ■ a%(x n ) (A™(y)A n p {z)) . (5.50) 



Nonpcrturbative calculational method in quantum held theory 



18 



If we suppose that k\ = 2 and fa = 4 then the rhs of equations l|5.30[) (|5.31|) will be multiplied on 2 and 
equations (|5.32|) (|5.34|) (|5.35|) will be multiplied on 4. Consequently the full Lagrangian will be 



- ^(0>)^)) 2 + ^a>V^ c 



(5.51) 



This Lagrangian we will investigate numerically on the next sections. 

Our final result given in l|5.49|l depends on several crucial assumptions (e.g. the existence of the negative 

2 

mass term — 12 ^-((f> a (p a )). In the next section we will investigate a non-topological flux tube solution to this 
effective Lagrangian. 

5.4 Flux tube equations 

We will start from the SU(2) Yang - Mills - Higgs field equations Ij5.49|l with broken gauge symmetry 

V v h afil/ = ge abc ct) b V^<l) c - (m 2 ) ab a b », (5.52) 
= -\r{4> b 4> b ~4L) (5.53) 

here h aflv = <9 M a" — B v a a + ge abc a h a^ is the field tensor for the SU(2) gauge potential a°; a,b,c = 1,2,3 are 
the color indices; D v \- ■ - ] a = d v [- ■ - ] a + ge abc a b [- ■ -} c is the gauge derivative; (f> a is the Higgs field; A, g and 
</>oo are some constants; (to 2 )" 6 is a masses matrix which destroys the gauge invariance of the Yang - Mills - 

Higgs theory, here we choose (m 2 ) ab = diag {to 2 , to 2 ., 0}. With a standard SU(2)-Higgs symmetry breaking 
Lagrangian the masses developed would be identical, and would be related to the gauge coupling and vacuum 
expectation of the field - m a ~ 4>oo9- In the present case we are postulating that the masses of the two gauge 
bosons are different. This was done in order to find electric flux tube solutions for the system (|5.52|l l|5.53|l . 
For equal masses electric flux tubes did not exist for the ansatz used. In numerically solving (|5.52|l l|5.53|l the 
two mass come out close to one another, but are not equal. This splitting of the gauge bosons mass over and 
above the Higgs mechanism explicitly breaks the gauge invariance, and at present we have no good justification 
for this step other than the fact that it leads to the existence of electric flux tubes of the kind that are very 
important in confinement. 

We will use the following ansatz 

V ^ 2/ \ V (P) ,3, \ HP) (r r A \ 

<w) = — ; <w) = — ; <p(p) = — (5-54) 

9 9 9 

here z, p, ip are cylindrical coordinates. Substituting this into the Yang - Mills - Higgs equations l|2.1|l (|5.53|) 

gives us 

/" + - = f(4> 2 + v 2 -m\), (5.55) 

" '' «(0 2 -/ 2 ~m 2 ), (5.56) 

4>" + ~ = 0[-/ 2 + « 2 + A(^ 2 -O] (5.57) 
x 

here we redefined cj>/4> -> <j>, f/<f> -> /, v/<f> -> v, (f>oo/(f>o -* <f>oo, m lt2 /(f>a -> mi, 2 , P<fio —> x; <f> = <j>(0). 
Similar equations with the presence of A 1 ^ and without masses mi ; 2 were investigated in ref's [H3] where it 
was shown that the corresponding solutions becomes singular either on a finite distance from the axis p = or 
on the infinity. The solutions with the masses m 12 and a magnetic field H z were obtained in ref. and the 
result is that the following versions of the flux tube exist: (1) the flux tube filled with electric/magnetic fields 
on the background of an external constant magnetic/electric field; (2) the Nielsen-Olesen flux tube dressed with 
transversal color electric and magnetic fields. 

5.5 Numerical algorithm 

As in section 0] we will solve the equations set I|5.55ll - (l5.57ll with the iterative procedure which is described in 
more detail in ref. [HZj- Eq. H5.56(l for the i th step has the following form 



x 



fi-i - mi A , (5.58) 
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here the functions and were defined at the (i — l) th step and the null approximation for the functions 
/(x) and (j){x) is o (x) = 1.3-0.3/cosh 2 (x/4) and / (x) = 0.2/ cosh 2 (x). A gain these forms are chosen because 
they are already similar to the form of solution which is being sought. The numerical investigation for eq. (|5.58|l 
shows that there is a number m\ $ for which: for m,2,i > m?, j the solution is singular Vi(x) — > — oo by x — > oo 
and for TO2,i < to*, i the solution is also singular Vi(x) — > +00 by x — > 00. In Fig. 0the singular solutions are 
presented. The asymptotical form of the equation for the singular solution is 




Figure 7: The typical singular solution Vi(x) for equation (|5.58(l 



v' 

<« + -j- « " ™2,i) (5-59) 

and corresponding asymptotical singular solutions are 

e x V<t>i-i-™2,i 

« ±a± p= (5.60) 

\Jx 

where a± are some constants. This indicates that there may exist some value, m 2 .i = m 2^ f° r which the 
solution Vi{x) is regular. The number to*, i is found with the method of the iterative approximation as for m* 
and /j,* in the glueball case. It is necessary that the functions fi-i(x) — > and 4>i-i(x) — ► 0^ j_i as x — > 00 in 
order to have acceptable solutions. 

Eq. H5.55fl a t this step has the following form 

f" +~ =/i(^_i+«i -™y . (5-6i) 

here the function <^_i was defined on (z — l)" 1 step and the function Uj is the solution of eq. I|5.58fl . The 
numerical investigation for eq. Ij5.61|l shows that there is a number m\ i for which: for roi.j > m* i the solution 
is singular /j(x) — > —00 by x — ► 00 and for mi t i < to* i the solution is also singular /j(x) — * +00 by x — > 00. 
The asymptotic singular behavior for /, (x) is the same as in the previous case. This implies that there may be 
some particular value mi^ = to* i for which the solution fi(x) is regular. The number m\ i is found with the 
method of iterative approximation. It is necessary that the functions Vi(x) — > and <fo_i(x) — > (jf^ i as x — > 00 
in order to have acceptable solutions. 

The next step in the i th iteration is solving of eq. I|5.57|l 

4>'l + « 4>i [-fi +vf + \ (tf - c/tj] (5.62) 

here the functions fi and Vi are the solutions of eqns. (|5.58|) (|5.61l) . The numerical investigation for this 
equation shows that there is a number (i>* ooi for which: for 4>oo,i > (jf^ i the function cf>i(x) oscillates with 
decreasing amplitude, for <poo,i < ^to » the function 4>i{x) — > +00 by x — * xq. The singular solution in this case 
is similar to the solution presented on Fig. ^ The number is intermediate between the values that give the 
singular and oscillating solutions. For this value it is possible that a regular solution exists. To have acceptable 
behavior for the solution 4>i(x) it is necessary that the functions Wj(x) — > and fi(x) — > as x — > 00. 

For the initial values we chose the following values: /(0) = 0.2, u(0) = 0.5 and (j>(0) = 1.0, and /'(0) = 
u'(0) = 0'(O) = 0. The set of equations l|5.55|l - H5.57|) has one independent parameter A, which we set to A = 0.2 
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i 


1 


2 


3 


4 


m U 


1.194482.. . 


1.20214.. . 


1.20106... 


1.2011... 


m h 


1.1405457... 


1.1517068... 


1.151087... 


1.150935... 


Cm 


1.2360886... 


1.23189178.. . 


1.23253772... 


1.23260797.. . 



Tabic 2: The iterative parameters ^ i and 



to perform the numerical solution of the equations. The other parameters of the equations (ml 2 ^ i and </>^ i ) 
are determined in the process of solving the equations. 

The iterative process described above gives m% 2 -, i and (p^^ presented on Tabled The functions fi{x), Vi(x) 
and (j>i(x) for i — 1,2,3,4 are presented in figs. 15119"! and 1 101 As in the glueball analysis we have only numerically 
integrated the solution out to finite x so it is not clear if this solution is truely well behaved as x — > oo. For the 
range of x where we have found the solution it behaves as desired. 




Figure 8: The first four iterations for the A] gauge 
potential component. 




Figure 9: The first four iterations for the A\ gauge 
potential component. 
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Figure 10: The first three iterations for the 3 
scalar field. 
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Figure 11: The color electric E^(x), E^(x) and 
magnetic H^(x) fields 



In terms of the ansatz functions the color electric and magnetic fields are 



El(x) = F t 3 z = 
El(x) = Fl p 



, 3 _ f(x)v(x) 



9 

fix) 
9 



H 2 (x) = xe vpz F 2 "* = -x^ 



(5.63) 
(5.64) 
(5.65) 



These shape of these fields is presented in fig. ^2 
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From eqs. H5.55M5.57jl one can see that the asymptotic behavior of the regular solutions f*(x), v*(x) and 
cj>*(x) is 



fW = fo 
v*(x) 

4>*(x) 



v - 



ix 



60 7x^ + 

V X 



(5.66) 
(5.67) 
(5.68) 



where fo, vq and cj)o axe constants. 

From the Lagrangian in (|5.49|) the energy density e(x) in terms of the ansatz functions is given by 

2g 2 e(x) = f'\x) + v' 2 {x) + 4>' 2 {x) + v 2 {x)f(x) + v 2 {x)cf 2 (x) + 



f\x)^{x) + m*J\x) - m* 2 v\x) + - 



(5.69) 



Its shape is shown in Fig. 1121 The exponentially decreasing asymptotic forms indicate that the energy density 
should go to zero at the infinity, and thus the total energy of the system should also be finite. Since we have 
only been able to integrate the solution to finite x we can not say absolutely that the total energy is finite. 
In the figs. El an d El the profiles of the linear energy density (the integral of the energy density over the p 
coordinate) is shown for different A and (f>Q. The linear energy density is given by 



2irpe(p)dp 



(5.70) 




Figure 12: The energy density e(x) 




Figure 13: The profile e(A) with the different O - 



Figure 14: The profile e(<?!>o) with the different A. 
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The flux of the longitudinal electric field is 



<T> = / E 3 z ds = 2tt p f ^ V ^ dp = — (°° xf(x)v(x)dx < oo. 
Jo 9 9 Jo 



(5.71) 



In the figs. 1151181 the profiles of the electric flux for different A and <po is given. Again the electric flux was 
calculated only out to some finite x. Eq. I|5.71|l assumes that the ansatz functions remain well behaved (i.e. 
approach zero exponentially as x — * oo) for larger x. 
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Figure 15: The profile $(A) with the different 4>q. 



Figure 16: The profile ^(A) with the different vq. 




Figure 17: The profile $(</>o) with the different 4>q. Figure 18: The profile 3>(0o) with the different ^o- 

We see that the flux $ depends on the parameters A, v(0) and 0(0) since the solutions f(x) and v(x) depend 
on these parameters. This is in contrast with such topological solutions as 't Hooft-Polyakov monopole and 
Nielsen-Olesen flux tube which have fixed fluxes. This is the indication of the fact that the colored flux tube 
solution is not the consequence of a first order differential equation (the Bogomolnyi Equation |?U]) as is the 
case for the topological solutions. The colored flux tube is a dynamical object not topological one. This 
raises questions about its stability, since this solution is then no longer guranteed to be stable via topological 
arguments. 



6 Discussion and conclusions 

Now we summarize the nonperturbative quantization method of this paper and enumerate some features of 
nonperturbative quantum field theory: 

• The existence and character of stationary, extended objects, such as flux tubes or glueballs must be 
explained using nonperturbative QFT . Such extended objects can not be understood within perturbative 
QFT as they are not simply a cloud of quanta. These objects share some similarity to a stationary 
turbulent liquid: there are fluctuations in any point within the liquid and simultaneously there exist 
averaged quantites (velocity, for example) which are correlated at any point as a consequence of the 
equations of motion for the turbulent liquid. 
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• It is well known that the ordinary commutations relations like 

n(x),<j)(y) = i5{x-y) 



(6.1) 



are correct only for linear field theory. The algebra of strongly interacting fields are more complicated and 
unknown. The algebra of nonperturbative field operators may be less singular then the algebra of linear 
fields in the sense that the Green's functions for the nonpertubative field operators may be more smooth 
and well behaved. As a consequence the fields at different spatial points but at one time will be correlated 
in contrast with the second (linear) case. The physical explanation for this is that in the linear case the 
quantized field at different points can be correlated only by exchanging quanta which move with some 
speed, but in the nonperturbative case these fields are correlated through nonlinear equations (similar to 
the correlation in a turbulent liquid) . 

• Nonperturbative QFT may be important for string theory. String theory says nothing about any possible 
inner structure of the string. By studying flux tube solutions such as those presented here one may obtain 
some insight at making a guess as to possible inner structure for strings. Here we would like to cite 
Polyakov |71| : "... It might be possible to visualize the superstrings as flux lines of some unknown gauge 
theory ..." 

As a final remark one can think to apply the present nonperturbative method to other strongly coupled, 
nonlinear systems such as gravity. In this way one could see if some insight into quantum gravity might not be 
gleaned. 
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A The flux tube in Euclidean spacetime 

In the Euclidean spacetime the equations (|5.55(l - (|5.57(l can be simplified. Making the transition to Euclidean 
spacetime via f(x) — > if(x) we have the following equations 



/' 



x 



v H 

x 

4> ^ — 

X 



f (0 2 + v 2 - ml) , 
v {4? + f- m 2 2 ) , 
4>[f + v 2 + \(4> 2 - ( g )] 



(A.l) 
(A.2) 
(A.3) 



For simplicity one can consider the case f(x) — v(x) and mi = — m which then leads to only two equations 



f"+ J - = f(4> 2 + f-m 2 ), 
<!>" + - = 0[2/ 2 + A(0 2 -<O 



(A.4) 
(A.5) 



The numerical method for obtaining a solution is similar to the methods applied in sections [S] and the results 
are presented in figs. ^5] , |3U] and Table [3] In fig [21 the profile of the color electric and magnetic fields are 



i 


i 


2 


3 


4 


TO* 


0.9436768025... 


0.79033... 


0.81214284.. . 


0.8139867... 


Cm 


0.9175.. . 


1.2248.. . 


1.1186... 


1.12185.. . 



Table 3: The iterative parameters to* and (jf^ i . 

presented. In figs. 1221231 we show the flux $ of color electric field 

$ = / Elds - 2t: r p fipHf>) dp = — r xf(x)v(x)dx. 
J Jo 9 9 Jo 

through the tube with respect to parameters A and <pQ. 



(A.6) 
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Figure 19: The iterative solutions /, (x) for the flux Figure 20: The iterative solutions 4>i(x) for the 

tube in Euclidean spacetime. flux tube in Euclidean spacetime. 



0.4- 




Figure 21: The profiles of the color electric and magnetic fields. 



B The numerical calculations of the ID A0 4 soliton 

As a test of the numerical method of solving the nonlinear equations l|4.33|l i|4.34[) we will apply this method to 
a known solution namely the A</> soliton in ID. The equation for this system is 

g=," = y (l- y 2 ) , (B.l) 
and the solution is _ 

•w - =E- < a2) 

We rewrite l|B.l|l in the form of the Schrodinger equation 

-y" + yV eff = -Ay (B.3) 
where V e ff — —y 2 and A = 1. We will solve this equation by an iterative procedure. We start with the equation 

-y" + vi{-vl) = -Am (B.4) 

where y\ (x) and Ai are the first approximations of the solutions. For the numerical solution we choose the null 
approximation as 

Vo = ^ r . (B.5) 
cosh(f) 

A typical solution for arbitrary values of the parameter Ai is presented in fig l24l This picture shows that 
there is a value AJ for which the solution is exceptional one. One can find this exceptional solution choosing 
the appropriate value of the "energy level" A* . After which an exceptional solution y\ (x) is substituted into 
equation for the second approximation y2(x) 

-y'i - V2 {vlf = -A22/2 (B.6) 
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Figure 22: The flux of the color electric field with 
respect to A with different 4>q. 
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Figure 24: The singular solutions for the soliton 
equation. 
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Figure 23: The flux of the color electric field with 
respect to 0o and with different A. 




Figure 25: The iterative functions 2/1,2, 3- 



and so on. The result is presented in Tableland figEl One can see that A* — > 1 and y*(x) is convergent to 



i 


1 


2 


3 


A* 


1.406929666... 


1.148564915.. . 


1.03968278... 



Table 4: The iterative values of the parameter A. 
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